Some geometric critical exponents for percolation and the random-cluster model 
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We introduce several infinite families of new critical exponents for the random-cluster model 
and present scaling arguments relating them to the fe-arm exponents. We then present Monte 
Carlo simulations confirming these predictions. These new exponents provide a convenient way to 
determine fc-arm exponents from Monte Carlo simulations. An understanding of these exponents also 
leads to a radically improved implementation of the Sweeny Monte Carlo algorithm. In addition, our 
Monte Carlo data allow us to conjecture an exact expression for the shortest-path fractal dimension 

d m in in two dimensions: d m ln = (g + 2)(g + 18)/(32p) where g is the Coulomb-gas coupling, related 
to the cluster fugacity q via q — 2 + 2cos(g7r/2) with 2 < g < 4. 

PACS numbers: 05.50.+q, 05.10.Ln, 64.60.ah, 64.60.De 



The random-cluster model [l| is a correlated bond- 
percolation model that plays a central role in the the- 
ory of critical phenomena, especially in two dimensions 
where it arises in recent developments of conformal field 
theory [2j via its connection with Schramm-Loewner evo- 
lution (SLE) 0,0]. To each bond configuration A C E 
of a given graph G = (V,E), the random-cluster model 
assigns a weight proportional to q k ( A ">v\ A \ where k{A) is 
the number of clusters (connected components), \A\ is 
the number of bonds, and v,q > are parameters (see 
Fig. [T]). For q — 1 the random-cluster model reduces to 
independent bond percolation [j| , while for integer q > 2 
it provides a graphical representation of the g-state fer- 
romagnetic Potts model |6]. Furthermore, the q — > 
limit corresponds to uniform spanning trees, provided 
that v oc q a with < a < 1; this applies in particu- 
lar to the critical square-lattice random-cluster model, 
for which v = Jq [7]. The random-cluster model thus 
provides an extension of both percolation and the Potts 
model that allows all positive values of q, integer or non- 
integer, to be studied within a unified framework. 

When the graph G is planar, it is useful to map the 
bond configurations to loop configurations 0], as illus- 
trated in Fig. [T] The loop configurations are drawn on 
the medial graph, the vertices of which correspond to 
the edges of the original graph. The medial graph of 
the square lattice is again a square lattice, rotated 45°. 
Each unoccupied edge of the original lattice is crossed by 
precisely two loop arcs, while occupied edges are crossed 
by none. The continuum limits of such loops are of cen- 
tral interest in studies of SLE The outermost loop 
bounding a cluster defines its hull Q. 

In this Rapid Communication we shall define and 




Figure 1: Typical bond configuration A C E on a finite sub- 
graph of the square lattice: here \A\ = 31 and k(A) — 4. 
Also shown are the corresponding medial-lattice loops (inte- 
riors shaded for clarity). Here C m i n , 2 = |CxJ = 11; choosing 



ei = Xix 2 and e 2 = x 2 x 3 yields M minj2 = \Mx 2 ,e 2 1 = 16; and 
Tmin,2 = T X1 = 4. The graph-theoretic distance (shortest- 
path length) between x 2 and X3 is 3. 



study several infinite families of critical exponents for the 
random-cluster model, related to cluster size, hull length 
and shortest path. Some of these exponents have been 
studied previously for spanning trees or percolation only, 
while others appear to be entirely new. We shall give 
scaling arguments determining all but one of these ex- 
ponents as a function of q in the two-dimensional (2D) 
case, and Monte Carlo simulations confirming these pre- 
dictions. The remaining undetermined exponent is the 
shortest-path fractal dimension d m ; n .10], which relates 
the shortest-path length and Euclidean distance between 
two vertices on a cluster. However, our Monte Carlo data 
lead us to conjecture the exact formula 



= (g + 2)(g + W)/(32g) 



(1) 



2 



where q = 2 + 2 cos(gw/2) and g £ [2,4] is the Coulomb- 



gas coupling To our knowledge, d m i n has not previ- 



ously been studied for q ^ 0, 1. 

Our original motivation for studying these exponents 
arose out of practical concerns related to Monte Carlo al- 
gorithms, as we discuss below. However, we subsequently 
discovered a relationship between them and the fc-arm 
exponents Xk @, which have proved to be of fundamen- 
tal importance in critical phenomena, not least within 
the context of rigorous studies of percolation [l2| . The 
exponents we study here are defined via the scaling of 
very natural graphical observables, which are easily mea- 
sured in Monte Carlo simulations. Therefore, the rela- 
tionship between these exponents and the /c-arm expo- 
nents provides a new and convenient way to numerically 
estimate Xk, which can be used also in three dimensions 
(3D) where no exact expressions for Xk are known. This 
has allowed us, for example, to numerically verify in 3D 
that Xi equals the thermal dimension Xt when q = 1, a 
result previously tested only in 2D. 

Motivation. Monte Carlo simulations are an essen- 
tial tool in statistical mechanics, but they typically suffer 
from critical slowing- down [l3j ]: the autocorrelation (re- 
laxation) time t diverges as a critical point is approached, 
most often as a power law r ~ £ z , where £ is the spatial 
correlation length and z is a dynamic critical exponent. 

The Sweeny algorithm 114] is a local single-bond up- 
date dynamics for the random-cluster model; for < 
g < 1 it is the only known general algorithm for this 
model. We have recently shown [lf| that the Sweeny al- 
gorithm has unusually weak critical slowing-down, fairly 
close to the theoretical lower bound z > ajv JjJ. Fur- 
thermore, it exhibits (especially for small q) the surpris- 
ing phenomenon of "critical speeding-up" [15| , in which 
suitable global observables exhibit significant decorrela- 
tion on time scales much less than one sweep (namely, 
L w hits for some w < d); this makes the algorithm po- 
tentially very efficient. 

The main obstacle to the use of the Sweeny algorithm 
is the need (when g ^ 1) to determine how updating a 
bond xy affects k(A): this potentially requires travers- 
ing an entire cluster to determine whether x and y are 
connected. Given a lattice site x, we write C x for the 
cluster containing x and |C X | for the number of sites in 
it. Since |C X | has mean ~ L 1 ^ near the critical point, 
the connectivity check threatens to impose a "computa- 
tional critical slowing-down" that would more than out- 
weigh the good "physical" behavior of the Sweeny dy- 
namics. There do exist sophisticated algorithms in the 
computer-science literature fl7j j for performing such con- 
nectivity checking dynamically, which have been proven 
to be (asymptotically) very efficient, but their complexity 
appears prohibitive for use in practical simulations. 

Our interest in the present project began with a 
simple idea for reducing this computational slowing- 
down without the need for complex data structures 



or algorithms: namely, perform simultaneous breadth- 
first searches starting at both endpoints x and y, and 
stop when one of the clusters has been fully visited 
or the clusters merge. In the first case this takes a 
time min(|C x |, |C y |), and in the second case a time B s 
(the number of sites visited in breadth-first search un- 
til merger). A natural question is therefore to the de- 
termine the critical behavior of C m i n! 2 and B s , where 
Cmin,2 = rnin(|C x |, |C y |) if C x and C y are distinct and 
otherwise. We will provide here a scaling argument sug- 
gesting that both scale (in mean) as L dp ~ X2 <C L 1 ^ , 
where dp = d — j3/v is the cluster fractal dimension and 
X2 is the 2-arm exponent. We will then verify this pre- 
diction numerically. 

In two dimensions an even more efficient procedure 
is to simultaneously follow the medial-lattice loops sur- 
rounding x and y. If the two loops are distinct, this 
takes a time min(|.M XiXy |, |-A/f yjXy |), where 7W x , e is the 
loop on the medial lattice that winds around x through 
the edge e if e is unoccupied, and M.^ e — other- 
wise (see Fig. [I]). This naturally leads to the question 
of the scaling of M m j ni 2, defined analogously to C m - ln .2 
as M minj2 = rnin(|.M x , xy |, |-M y , xy |) if -M x , xy and -M y , xy 
are distinct and otherwise. We will provide a scaling 
argument, and confirm numerically, that -M m j nj 2 scales 
(in mean) as ~ L d H-x 2 wnere g] H = \ _|_ 2/g is the 
hull fractal dimension 0, [3]. It follows that compu- 
tational critical slowing-down is completely absent for 
q > 4 cos 2 (tt^/2/3) rs 2.811520. 

Definition of exponents. These scaling results for 
Cmin,2 and M m i n> 2 can be generalized in a very natu- 
ral way. We shall consider a variety of positive-integer- 
valued observables O; for each one we expect that its 
probability distribution obeys a scaling law ¥(0 = s) ~ 
s~vo ( w ith ipQ > 1) for large s at criticality in infinite 
volume, or more generally ~ s~^° Fo(s/£ d ° , s/L d °) near 
criticality in large finite volume, where Fq is a scaling 
function. Our goal is to determine, for each O, the decay 
exponent ipo and the fractal dimension do . Note that at 
criticality in finite volume, (O n ) ~ L(n+i-^ )d + const 
+ corrections to scaling as L — > oo (or ~ log'L if 
n + 1 — tpo = 0); we shall use this fact in our Monte 
Carlo determinations of ipo and do- 

Fix nearby sites Xi,...,Xfe, and let C m i rii fc = 
min(|C Xl |, . . . , |C X J) if these clusters are all distinct, and 
otherwise. We expect that all these observables have 
fractal dimension dc mln k equal to the cluster fractal di- 
mension dp = d — /3/v. Moreover, standard hyperscaling 
arguments [H[ give ^c min i = d/dF (the usual notation 
is t = "0c min i + !)■ Our goal is to determine V , c min fc 
for k > 2. To our knowledge these exponents are new, 
except V , c min 2 — if f° r spanning trees [l9[ . 

For two-dimensional lattices, choose for each site 
Xi a bond incident on it, and let Af m i nj /c ~ 
min(|A / l Xl!ei |, . . . , \M Xk ^ k |) if these loops are all distinct, 
and otherwise. We expect that all these observables 



3 



have fractal dimension c?A/ min k equal to the hull fractal 
dimension dn ', and standard hyperscaling arguments give 
V'A/min i = d/dfi- We aim to determine ^M min k f° r k > 2. 

Now let T x denote the maximum graph-theoretic dis- 
tance from x to any site in C x , and define Tmin.fc = 
min(T Xl , . . . , T Xfc ) if the clusters C Xi are all distinct, and 
otherwise (sec Fig. [I}. We expect that all these observ- 
ables have fractal dimension dx alin k equal to the shortest- 
path fractal dimension d m i n . We aim to determine V'T min k 
for k > 1, as well as d m i n . 

Finally, consider a pair of sites separated by a distance 
~ L, say x = and x = aL where a 6 M. d , and let S a L 
be the length of the shortest path connecting these sites if 
one exists, and otherwise. Since P(0 f-> aL) ~ L~ 2 P/ V , 
we expect that ({S aL ) n ) ~ L nd ™- 20 /". We will use 
((S a L) n ) to numerically determine d m i n in 2D, leading 
to the conjecture (p}. 

Scaling arguments. Let Pk(R) be the probability that, 
in an annulus of inner radius r ~ 0(1) and outer radius 
R, the inner circle is connected to the outer one by at 
least k distinct clusters. The fc-arm exponent Xk char- 
acterizes the large- R asymptotics of pk (R) at criticality: 
Pk(R) ~ R~ Xk as R — V oo. In particular, x\ = fl/v. In 
two dimensions it is known 0, [III [2(| that 



xx = (g-2)(6-g)/(8g) 

(. 9 /8)fc 2 -( 5 -4) 2 /(8<?) 



for k > 2 



(2a) 
(2b) 



Now fix nearby sites Xi,. . . ,x&, and let Cfc(s) be the 
probability that xi, . . . ,x* belong to k distinct clusters, 
each of which contains at least s sites. The correspon- 
dence s ~ R d " suggests that Ck{s) ~ pk{s 1 ^ dF ). Since 
Cfe(s) = P(C mirl ,fc > s), we predict if>c min , k = x k /d F + 1 
and hence ((C m in,fc)") ~ L ndF ~ Xk + const. This agrees 
with standard hyperscaling for k = 1; we will test it 
for k > 2. A similar argument suggests that (B™) ~ 

Analogous reasoning predicts V>Af min k = £fc / <ifr + 1 and 
V'Tmin fc = Xk/ drain + 1- For q = 1 this argument for 



shortest-path scaling is due to Ziff 2l|. In 2D it is known 
0, [53 that d H = 1 + 2/g, but to our knowledge (J) is 
the only known conjecture for d min , even for q = 1. 

Monte Carlo simulations. We simulated the two- 
dimensional random-cluster model at criticality for q = 
0, 0.01, 0.25, 0.5, 1, 1.5, 2, 3, 3.5 ^ and 4 < L < 1024 (pe- 




Figure 2: (Color online) Finite-size-scaling plot showing 
s X2/dF P(C m in,2 > s) versus s/L dF for g = 1 and 64 < L < 
4096. 



(the smallest L value included in the fit) until the % 2 was 
reasonable. The error bar is a subjective 68% confidence 
limit that includes both statistical error and systematic 
error due to unincluded corrections to scaling. 

Our results for (C m i n ,2) and (Af m i n ^) are presented 
in Table Q] complete results for (C" in 2 




and 

(M™ in 2 ) with n = 1,2,3 can be found in [26j. The agree- 
ment with the predicted exponents is excellent, except 
where the exponent is very negative and hence possibly 
overshadowed by correction-to-scaling terms. A finitc- 
size-scaling plot for C m i n! 2 and q = 1 is shown in Fig. [5] 
and exhibits excellent collapse. 

Next we studied S a L for a = (i, 0) in order to estimate 
d m in- see Table U and 26(. Our result for q — 1 is compat- 
ible with Grassberger's [13] estimate d m ; n = 1.1306(3). 
The corrections to scaling are very strong for these ob- 
servables, and our error bars are dominated by our as- 
sessment of the likely systematic error from such correc- 
tions. It would be very useful (but also very expensive) 
to obtain data at larger values of L. 

Our results for d m j n are consistent with the simple 
formula (HJ: see Table U and Fig. [3] This formula has 
the nice property that d m ; n is monotone decreasing for 
2 < g < 6 and reaches d m in = 1 precisely at g = 6, in 
accordance with the idea [28[ that clusters become more 
compact as g grows. It also agrees with the known fact 



riodic boundary conditions) by the Sweeny algorithm H that d min = | at q = . Indeed, if one seeks a formula 



when < q < 1 and the Chayes-Machta algorithm 23] 
when q > 1. For 0.25 < q < 2 we also have data at 
L = 2048. For q — 1 we used cluster-growth algorithms 
to handle 4 < L < 4096. For q = we used Wilson's al- 
gorithm [24| to generate spanning trees from loop-erased 
random walk [25j. The total CPU time used in these 
simulations was approximately 66 years using a 3.2 GHz 
Xeon EM64T processor. 

For each observable O, we fit (O n ) for n = 1,2,3 to 
the Ansatze aL p , aL p + bL p ~ A and aL p + c, varying L m i n 



of the "Coulomb-gas" form <i m i n = F(g) = Ag + B + C/g 
11] and imposes the constraints F(2) = |, F(6) = 1, 
F'(6) = 0, then the unique solution is (p}. 

Our results for Xim^i and T m i nj 2 are similar to those 
shown in Table U and will be reported elsewhere. 

We are grateful to Gary Reich for a discussion that in- 
spired this work; to Wenan Guo for participation in early 
stages of this project; to Russ Lyons, Robin Pemantle and 
David Wilson for helpful correspondence; and to NYU 
ITS for use of their computer cluster. This research was 
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liiiat"i^"i'f~'\7" 
Vc^J LtttllLl L y 


q= 0.01 0.25 0.5 1 1.5 2 3 3.5 


(C m in,2) 


in ( T\ 1 1 TYl 1 
jj I I1U.II1 . j 

p (pred.) 


1.25000 1.18007 0.92707 0.80768 0.64583 0.52298 0.41667 0.21667 0.10376 


(M min , 2 > 


p (num.) 
p (pred.) 


1.25002(4) 1.1639(5) 0.8527(5) 0.7036(3) 0.4994(7) 0.3444(4) 0.2087(2) -0.0525(25) -0.195(3) 
1.25000 1.16439 0.85222 0.70341 0.50000 0.34420 0.20833 -0.05000 -0.19748 


Sl/2 
2<9<4 


a!min (num.) 
a!min (conj.) 


1.24999(3) 1.2371(10) 1.1825(3) 1.1596(4) 1.1303(8) 1.1112(7) 1.0955(10) 1.0677(40) 1.0560(30) 
1.25000 1.23463 1.18211 1.15918 1.13021 1.10997 1.09375 1.06667 1.05343 



Table I: Numerical estimates versus theoretical predictions for exponents associated to {C m i n ,2}, (Af m i n ,2) and Sl/2- 



drain 


f predicted 




t numerical j 


1.2 






N^^critical (2 < g < 4) 


1.1 




1.0 


tricritical (6 > g > 4) 





1 2 3 4 



Figure 3: (Color online) Numerical estimates for d m i n to- 
gether with the conjectured exact formula d m m = (g + 2)(g + 
18)/(32g). 
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Detailed results for (C™ in 2 ) and (C™ in 3 ) with n = 
1,2,3 are presented here in Table I. The agreement with 
the predicted exponents is excellent, except where the 
exponent is very negative and hence possibly overshad- 
owed by correction-to-scaling terms. (Exponents p < 
always come from a fit to aL p + c.) For (C m i n ,3) at q = 0, 
where the predicted exponent is zero, we observe a clear 
alogL + b behavior: see Fig. 1. 

Detailed results for (M™ in 2 ) with n — 1,2,3 are re- 
ported in the top half of Table II. The predicted val- 
ues are again confirmed. We also studied M m i n ^ for 
the "dense" Fortuin-Kasteleyn (FK) clusters in which 
all bonds between sites in the same cluster are inserted. 
This modification leaves connectivities (hence Xk and dp) 
unchanged but is conjectured [1] to change the hull and 
shortest-path dimensions to those of FK clusters with 
coupling g' = 16/g. The results for M m ; ni2 are shown 
in the bottom half of Table II: they are consistent with 
our prediction over the entire range 4 < g' < 8. We re- 
mark that 4 < g' < 6 corresponds to the tricritical Potts 
model for 4 > q > 0, but no interpretation of 6 < g' < 8 
in terms of a spin model with positive weights seems to 
be known [2]. 



Results for (S™ L ) with a = (^,0) and n — 1,2,3 are 
presented in Table III. These data are consistent with the 
scaling (S™ L ) ~ ind m _in-W/ v an( j with the conjectured 

exact formula d m i n = (g + 2)(g + 18)/(32g). 



[1] Y. Deng, unpublished. See also Y. Deng, H.W.J. Blote and 
B. Nienhuis, Phys. Rev. E 69, 026123 (2004) for a related 
construction. 

[2] The critical 0(n) loop model with 2 > n > is known 
to correspond to the tricritical Potts model with q = n 2 
and hence 4 < g < 6 [B. Nienhuis, Physica A 177, 109 
(1991)]. The critical 0(n) loop model with > n > -2 
may correspond to 6 < g < 8, but it has negative weights. 




Quantity 


q= 0.01 0.25 0.5 1 1.5 2 3 3.5 


C*min,2 


p (num.) 
p (prcd.) 


1.25002(4) 1.1800(1) 0.9273(2) 0.8071(2) 0.6458(1) 0.5226(2) 0.4163(3) 0.214(2) 0.108(3) 
1.25000 1.18007 0.92707 0.80768 0.64583 0.52298 0.41667 0.21667 0.10376 


(C*min,2) 2 


p (num.) 
p (prcd.) 


3.25002(4) 3.1647(2) 2.8633(1) 2.7252(2) 2.5418(2) 2.4065(2) 2.2923(3) 2.0832(7) 1.9695(8) 
3.25000 3.16489 2.86336 2.72492 2.54167 2.40621 2.29167 2.08333 1.96999 


(C min , 2 ) 3 


p (num.) 
p (prcd.) 


5.24998(5) 5.1496(2) 4.7997(2) 4.6424(2) 4.4377(2) 4.2896(3) 4.1668(6) 3.950(2) 3.838(2) 
5.25000 5.14970 4.79965 4.64216 4.43750 4.28943 4.16667 3.95000 3.83621 


^min,3 


p (num.) 
p (prcd.) 


a\ogL + b -0.109(2) -0.532(9) -0.723(8) -0.961(8) -1.04(3) -0.98(2) -0.7(1) -0.6(1) 
0.00000 -0.10974 -0.52401 -0.72989 -1.02083 -1.25148 -1.45833 -1.86667 -2.10868 


(Cmin,3) 2 


p (num.) 
p (prcd.) 


2.0009(5) 1.8755(12) 1.415(3) 1.184(4) 0.867(5) 0.623(5) 0.38(3) 0.01(2) -0.29(9) 
2.00000 1.87508 1.41229 1.18735 0.87500 0.63174 0.41667 0.00000 -0.24245 


(C min ,3) 3 


p (num.) 
p (prcd.) 


4.0012(9) 3.860(2) 3.352(4) 3.100(6) 2.739(18) 2.494(11) 2.33(7) 1.73(7) 1.56(5) 
4.00000 3.85990 3.34857 3.10460 2.77083 2.51497 2.29167 1.86667 1.62378 



Table I: Numerical estimates versus theoretical predictions for exponents associated to C m i n ,fc with k = 2, 3. 



Quantity 


<? = 





0.01 


0.25 


0.5 1 


1.5 


2 


3 


3.5 


Mmin,2 


p (num.) 


1.25002(4) 1.1639(5) 


0.8527(5) 


0.7036(3) 0.4994(7) 


0.3444(4) 


0.2087(2) 


-0.0525(25) 


-0.195(3) 




p (pred.) 


1.25000 


1.16439 


0.85222 


0.70341 0.50000 


0.34420 


0.20833 


-0.05000 


-0.19748 


(M min , 2 ) 2 


p (num.) 


3.25002(4) 3.1333(4) 


2.7141(5) 


2.5164(4) 2.2498(4) 


2.0488(2) 


1.8749(3) 


1.5485(15) 


1.372(2) 




p (pred.) 


3.25000 


3.13353 


2.71364 


2.51638 2.25000 


2.04864 


1.87500 


1.55000 


1.36751 


(M m ,„,2) 3 


p (num.) 


5.24998(5) 5.1023(4) 


4.5756(6) 


4.3293(5) 4.0000(2) 


3.7531(4) 


3.5415(3) 


3.1484(16) 


2.938(3) 


2<<y<4 


p (pred.) 


5.25000 


5.10267 


4.57507 


4.32935 4.00000 


3.75308 


3.54167 


3.15000 


2.93250 




p (num.) 


0.496(3) 


0.4532(6) 


0.2818(10) 0.1973(5) 0.085(3) 


-0.006(5) 


-0.099(5) 


-0.24(3) 


-0.46(15) 




p (pred.) 


0.50000 


0.45322 


0.28100 


0.19795 0.08333 


-0.00535 


-0.08333 


-0.23333 


-0.31998 


(ACin,2) 2 


p (num.) 


1.748(3) 


1.7116(6) 


1.5718(4) 


1.5050(8) 1.4165(13) 1.351(2) 


1.296(4) 


1.11(4) 


1.02(7) 




p (pred.) 


1.75000 


1.71118 


1.57121 


1.50546 1.41667 


1.34955 


1.29167 


1.18333 


1.12250 


(^in, 2 ) 3 


p (num.) 


2.995(4) 


2.9688(10) 2.8612(20) 2.8136(8) 2.749(4) 


2.709(10) 


2.667(2) 


2.59(4) 


2.55(4) 


4 < g' < 8 


p (pred.) 


3.00000 


2.96914 


2.86143 


2.81297 2.75000 


2.70444 


2.66667 


2.60000 


2.56499 



Table II: Numerical estimates versus theoretical predictions for exponents associated to M m i nj 2. Top half of table refers to 
ordinary FK clusters; bottom half of table referes to "dense" FK clusters in which all bonds between sites in the same cluster 
are inserted. 



Quantity 


q= 0.01 0.25 0.5 1 1.5 2 3 3.5 


Sl/2 
(Sl/2) 2 
(Sl/2) A 
2<S<4 


fimin (num.) 
rfmin (num.) 
d min (num.) 

fimin (conj.) 


1.24999(3) 1.2371(10) 1.1825(3) 1.1596(4) 1.1303(8) 1.1112(7) 1.0955(10) 1.0677(40) 1.0560(30) 
1.25000(3) 1.2364(10) 1.1826(3) 1.1596(4) 1.1304(10) 1.1108(6) 1.0950(9) 1.0668(30) 1.0538(25) 
1.24999(4) 1.2359(10) 1.1825(3) 1.1595(3) 1.1304(10) 1.1106(3) 1.0949(9) 1.0662(30) 1.0525(23) 
1.25000 1.23463 1.18211 1.15918 1.13021 1.10997 1.09375 1.06667 1.05343 



Table III: Numerical estimates for the shortest-path exponent rf m i n as deduced from ((S , i / 2 ) n ) for n = 1, 2, 3, versus conjectured 
exact formula. 



